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In this work we have combined the generalized augmented space method introduced by one of 
us with the recursion method of Haydock et.al (GASR), within the framework of the local density 
functional based linear muffin-tin orbitals basis (TB-LMTO). Using this we have studied the effect 
of short-range ordering and clustering on the density of states, optical conductivity and reflectivity 
of 50-50 CuZn alloys. Our results are in good agreement with alternative techniques. We argue that 
the TB-LMTO-GASR is a feasible, efficient and quantitatively accurate computational technique 
for the study of environmental effects in disordered binary alloys. 
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I. INTRODUCTION 

Binary alloys involving equal proportions of a noble metal 
Cu and a divalent metal Zn have a stable low temperature 
/3-phase which sits immediately to the right of the pure 
face-centered cubic Cu phase in the alloy phase diagram 
0,0- This phase, called j3 brass, has a body-centered 
cubic structure. At high temperatures the alloy forms 
a disordered body-centered cubic structure. At around 
730K it orders into the B2 structure with two atoms per 
unit cell. The alloy satisfies the Hume-Rothery rules Q 
and has the same ratio of valence electrons to atoms. 
Jona and Marcus Q have shown from a density func- 
tional theory (DFT) based approach that within the lo- 
cal density approximation (LDA), it is the body-centered 
based B2 which is the stable ground state. They also 
showed that if we include the gradient corrections (GGA) 
then we get a tetragonal ground state lower in energy by 
0.1 mRy/atom. This is in contradiction with the latest 
experimental data. The alloying of face-centered cubic 
Cu with an equal amount of Zn leads to a body-centered 
stable phase. Zn has only one more electron than Cu. 
This is an interesting phenomenon. CuZn alloys also 
have anomalously high elastic anisotropy. This makes 
the theoretical study of CuZn an interesting exercise for 
a proposed theoretical technique. 

One of the earliest first-principles density functional 
based study of the electronic properties of CuZn was by 
Bansil Q. The authors had studied the complex bands 
of a-phase of CuZn using the Korringa-Kohn-Rostockcr 
(KKR) method coupled with the coherent potential ap- 
proximation (CPA) to take care of disorder. They com- 
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mented on the effects of charge transfer and lattice con- 
stants on the electronic structure. They found the elec- 
tronic distribution of this alloy to be of a split band kind 
with the centers of the Cu and Zn d-bands well sepa- 
rated from each other. Their Zn d-bands showed hardly 
any dispersion and were shown only schematically in their 
figures. In a later work Rowlands |5j generalized the CPA 
to a non-local version (NL-CPA) and studied the effects 
of short-range ordering in CuZn. Their technique was 
based on an idea of renormalization in reciprocal space 
suggested by Jarrell and Krishnamurthy Q . 
The order-disorder transition in CuZn is a classical ex- 
ample of a true second order transition. Several very 
early investigations on this alloy have been reviewed by 
Nix and Shockley Q and Guttman ||. These investi- 
gations, of course, were rather crude, since sophisticated 
approaches to deal with disordered alloys had really not 
been developed at that time. However, it was recognized 
that a knowledge of the short-range order correlations 
above the critical temperature should be of considerable 
interest. Early neutron scattering experiments were car- 
ried out on /?-brass by Walker and Keating 0. The 
Warren-Cowley short-range order parameter, defined by 
a(R) = 1 — Pab{R)/x, where x was the concentration of 
A and Pab (R) was the probability of finding an A atom 
at a distance of R from a B atom, was directly obtained 
from the diffuse scattering cross-section : 

^ = x{l-x){b A -b B ) 2 ^2a(R)f(K)exp(iK-R) 

R 

bAj^B were the scattering lengths of A and B atoms, and 
f(K) = exp(— C|i^| 2 ) was the attenuation factor arising 
from thermal vibrations and static strains. The exper- 
imental data for the short-range order parameter as a 
function of temperature are thus available to us. A Ising- 
like model using pair interactions was studied by Walker 
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and Chipman [lO( and the short-range order was theo- 
retically obtained. However, the pair interactions were 
simply fitted to the experimental values of the transition 
temperature Tq and in that sense it was an empirical 
theory. The experimental estimate of the nearest neigh- 
bour Warren-Cowley parameter was found to be varying 
between -0.171 to -0.182 at around 750K. 

In a later work using the much more sophisticated locally 
self-consistent Green function (LSGF) approach based on 
the tight-binding linear muffin-tin orbital (TB-LMTO) 
technique Abrikosov et.al. 0] studied CuZn alloys. The 
authors argued that earlier studies of the mixing en- 
thalpies of CuZn using the standard coherent poten- 
tial approximation approaches [l2^-|lfij showed signifi- 
cant discrepancies with experiment. The discrepancies 
were assumed to partly arise from the neglect of charge 
transfer effects and partly because of short ranged order- 
ing (SRO). The main thrust of this technique, which was 
based on an earlier idea of a locally self-consistent mul- 
tiple scattering (LSMS) by Wang et.al. 0], was to go 
beyond the CPA and include the effects of the immedi- 
ate environment of an atom in the solid. The LSMS gave 
an excellent theoretical estimate of the ordering energy 
in CuZn : 3.37 mRy/atom as compared to the exper- 
imental value of 3.5 mRy/atom. The LSGF approach 
correctly predicted ordering tendency in CuZn on lower- 
ing temperature and combining with a Cluster Variation- 
Connolly Williams (CVM-CW) obtained a value of the 
nearest neighbour Warren-Cowley SRO parameter a = 
-0.15. Subsequently Bruno et.al 0] proposed a modi- 
fication of the CPA including the local field effects and 
showed that charge transfer effects can be taken into ac- 
count as accurately as the O(N) methods just described. 
They applied their approach to the CuZn alloys. 

One of the earlier works on the optical property of CuZn 
alloy was the determination of the temperature variation 
of optical reflectivity by Muldawer |l9j. The author at- 
tempted to explain the color of the disordered <3-brass 
CuZn alloy via the internal photoelectric effect |20j . Al- 
though, the experimental data also contained the contri- 
bution from plasma oscillations, the author claimed that 
the optical reflectivity helps to explain the band picture 
of the alloys as a function of the inter-atomic spacin g. In 
order to explain the optical properties, Amar et.al. |2l|- 
[23l | studied the band structure of CuZn using the KKR 
method. However, they had used the virtual crystal ap- 
proximation, replacing the random potential seen by the 
electrons by an averaged one. This is now known to be 
particularly inaccurate for split band alloys. 

The above discussion was necessary to bring into fo- 
cus the following points : in the study of alloys like 
CuZn it would be interesting to address the effects of 
charge transfer and short-range ordering. In this com- 
munication we shall address exactly these two points. 
We shall propose the use of the augmented space recur- 



sion (ASR) coupled with the tight-binding linear muffin- 
tin orbitals basis (TB-LMTO) to study the effects 
of short-range ordering on both the electronic structure 
and the optical properties of /3-CuZn alloy at 50-50 com- 
position. We should like to stress here that the TB- 
LMTO-ASR addresses precisely these effects with accu- 
racy : the density functional self-consistent TB-LMTO 
takes care of the charge transfer, while the local environ- 
mental effects which are essential for the description of 
SRO are dealt with by the ASR. The TB-LMTO- ASR 
and its advantages has been extensively discussed earlier 
in a review by Mookerjee [2^ and in a scries of articles 
El ElimHIiil- We would like to refer the inter- 
ested readers to these for details. 



II. SPECTRAL FUNCTIONS, COMPLEX BANDS 
AND DENSITY OF STATES FOR 50-50 CUZN 

In this section we shall introduce the salient features of 
the ASR which will be required by us in our subsequent 
discussions. 

We shall start from a first principle TB-LMTO set of or- 
bitals 0, H3| in the most-localized representation. This 
is necessary, because the subsequent recursion requires 
a sparse representation of the Hamiltonian. The TB- 
LMTO second order tight-binding Hamiltonian H' 2 'is 
described by a set of potential parameters : Cr, E u r , 
A# and or which are characteristic of the atoms which 
sit on the lattice sites labelled by R, and a structure ma- 
trix Srri which is characteristic of the lattice on which 
the atoms sit. For a substitutionally disordered alloy, 
the structure matrix is not random but the potential pa- 
rameters are and can be described by a set of random 
occupation variables {ur}. We may write : 

Crl = Ct n R + Cf (1 - n R ) 

and similar expressions for the other potential parame- 
ters. The random site-occupation variables {ur} take 
values 1 and depending upon whether the muffin-tin 
labelled by R is occupied by A or B-type of atom. The 
atom sitting at {R} can either be of type A (ur = 1) 
with probability x or B (ur = 0) with probability y. 

In the absence of short-range order, the augmented space 
formalism associates with each random variable ur an 
operator M.r whose spectral density is its probability 
density. 

p(n R ) = — lim lm(U \ ({n R + iS)l - Mr)' 1 | Tr) 

7T (J — >0 

The operator M# acts on the "configuration space" of 
the variable ur, &r spanned by the configuration states 
I \r) and I [r). The augmented space theorem |2(| states 
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that a configuration average can be expressed as a matrix 
element in the "configuration space" of the disordered 
system : 

« A({n R }) » = < {0}|A({M«})|{0} > (1) 

where, 

A({M R }) = J ...J A({\ R })l[dP(X R ). 

P(Aij) is the spectral density of the self-adjoint opera- 
tor Mr, and the configuration state |{0}) is Hi I Tfl)- 
Applying (JTJ to the Green function we get : 

«G(k,z)» = (k®{0}|(zI-H^f 'ik® {0}). (2) 

where G and H^ 2 ) are operators which are matrices in 
angular momentum space, and the augmented k- space 
basis |k, L g) {0}) has the form 

(l/VAO J2 expHk • R)\R, L ® {0}). 

R 

The augmented space Hamiltonian H^ 2 ' is constructed 
from the TB-LMTO Hamiltonian H^ 2 ) by replacing each 



random variable n R by the operators M/j. It is an oper- 
ator in the augmented space \I> = 7i® JJ® &r. The ASF 
maps a disordered Hamiltonian described in a Hilbcrt 
space Tl onto an ordered Hamiltonian in an enlarged 
space ^P, where the space <J> is constructed as the outer 
product of the space H and configuration space $ of the 
random variables of the disordered Hamiltonian. The 
configuration space <& is of rank 2 N if there are N muffin- 
tin spheres in the system. Another way of looking at H*- 2 ' 
is to note that it is the collection of all possible Hamilto- 
nians for all possible configurations of the system. 
This equation is now exactly in the form in which recur- 
sion method may be applied. At this point we note that 
the above expression for the averaged G^^k, z) is exact. 

The recursion method addresses inversions of infinite ma- 
trices of the type associated with the Green function • 
Once a sparse representation of an operator in Hilbcrt 
space, H*- 2 - 1 , is known in a countable basis, the recursion 
method obtains an alternative basis in which the opera- 
tor becomes tridiagonal. This basis and the representa- 
tions of the operator in it are found recursively through 
a three-term recurrence relation. The spectral function 
is then obtained from the continued fraction : 



«G LL (k,z)» = @il = -ft (3) 

/?2i( k ) z-E L {k) -Ei(k,2) 



~ aiLik) 



z - a 2 i(k) 



z - ajvi(k) - r L (k, z) 



where T l (k, z) is the asymptotic part of the continued 
fraction. The approximation involved has to do with 
the termination of this continued fraction. The coeffi- 
cients are calculated exactly up to a finite number of 
steps {a n , (3 n } for n < N and the asymptotic part of the 
continued fraction is obtained from the initial set of co- 
efficients using the idea of Beer and Pettifor terminator 
[3^1 . Haydock and coworkers j3f| have carried out exten- 
sive studies of the errors involved and precise estimates 
are available in the literature. Haydock [36| has shown 
that if we carry out recursion exactly up to N steps, the 
resulting continued fraction maintains the first 2N mo- 
ments of the exact result. 



The self-energy Ei(k, z) arises because of scattering by 
the random potential fluctuations. 

The average spectral function <C Ak (E) ^S> is related to 
the averaged Green function in reciprocal space as : 

« A k (E) » = ]T « A kL (E) », 

L 

where 

< Awl(E) > = lim {Im < G LL (k 7 E - iS) >}. 

To obtain the complex bands for the alloy we fix a value 
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N 



FIG. 1: (left) Bands for pure Cu and Zn in bcc lattices with th 
line are for Cu and the full lines for Zn. (right) Complex bands 



for k and solve for 



£?i(k)-E L (k,£0 =0. 



The real part of the roots will give the position of the 
bands, while the imaginary part of roots will be propor- 
tional to the disorder induced broadening. Since the alloy 
is random, the bands always have finite lifetimes and are 
fuzzy. 

We have used this reciprocal space ASR to obtain the 
complex bands and spectral functions for the CuZn al- 
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FIG. 3: (bottom) Density of states of pure Zn (solid lines) and 
Cu (dashed lines)in the same bcc lattice as the 50-50 CuZn 
alloy, (centre) Density of states for ordered B2 50-50 CuZn 
alloy, (top) Density of states for the disordered bcc 50-50 
CuZn alloy. These results are comparable to the single-site 
CPA. 



FIG. 2: Spectral functions for the CuZn alloy for k-vectors 
along the r to iV direction in the Brillouin zone. 



